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Abstract 

Yukawa couplings between Higgs bosons and quarks in supersymmetric theories 
can be measured in the processes e + e~ — > QQ + Higgs. We have determined the 
cross sections of these processes in the minimal supersymmetric model including the 
complete set of next-to-leading order QCD corrections for all channels. 
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1. In the Standard Model (SM) and its supersymmetric extensions, the masses of elec- 
troweak gauge bosons, leptons, and quarks are generated by interactions with Higgs fields 
|Ij]. The Yukawa couplings between Higgs particles and fermions therefore grow with the 
masses Mf of the fermions. The couplings obey a universal scaling law = Mf/v in 
the SM, with v ~ 246 GeV being the ground-state value of the Higgs field. In supersym- 
metric models, which involve at least two Higgs doublets, the size of the Yukawa couplings 
is also set by the fermion masses, yet the relationship is more complex owing to the mixing 
among the Higgs fields. The Yukawa couplings of the two CP-even light/heavy Higgs 
bosons h/H and of the CP-odd Higgs boson A in the minimal supersymmetric extension 
of the Standard Model (MSSM) 0,f|, expressed in units of the SM couplings, are collected 
in Table T T . The parameter tan/3 = v%/vx is the ratio of the vacuum expectation values 
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Table 1: Higgs Yukawa couplings in the MSSM in units of the SM couplings Mf/v, and 
gauge couplings in units of the ZZH coupling Mz jv in the SM. 



of the two Higgs fields generating the masses of up- and down-type particles, and a is 
the mixing angle in the CP-even sector. In the decoupling limit, in which the light Higgs 
mass reaches the maximum value for a given parameter tan/3, the h Yukawa couplings 
approach the SM values. In general, they are suppressed for up-type fermions and en- 
hanced for down-type fermions; the enhancement increases with tan f3 and can therefore 
be very strong. 

Higgs radiation off top or off bottom quarks in e + e~ collisions, 

e + e" -> QQh/H/A [Q = t,b], (1) 

lends itself as a suitable process for measuring the Yukawa couplings in supersymmetric 
theories [BJ, particularly for the light Higgs boson h and for moderately heavy Higgs bosons 
H and A. In this letter we present the cross sections for these processes including the next- 
to-leading order (NLO) QCD corrections. This problem has been solved before in the SM, 
where the QCD corrections significantly increase the cross section at moderate collider 
energies |7|||. Moreover, it has been shown in detailed simulations that the SM top-quark 
Yukawa coupling can be measured with an accuracy close to 5% M (see also Ref. |TD[) 



in high-luminosity operations of e e linear colliders. For the MSSM, first steps towards 
an NLO analysis have been made in Ref. [ LT| , adopting either a global K factor from 



1 For large tan/3, SUSY loop corrections to the 66-Higgs vertices become very important ||. In 
the present QCD analysis, these corrections can be absorbed into effective Higgs Yukawa couplings, 
since emission and reabsorption of virtual heavy SUSY particles is confined to small space-time regions 
compared with QCD subprocesses involving massless gluons. 
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SM analyses or applying the QCD corrections to the cross sections for double resonance 
production in the narrow-width approximation. In the present analysis we have calculated 
the complete set of 0(a s ) QCD corrections to all subchannels in (|l|) systematically. The 
large number of interfering subchannels in supersymmetric theories renders this program 
more complex than the corresponding calculation in the SM, cf. Fig. |l|, in particular since 
the relative weight of the subchannels varies over the supersymmetric parameter space 
and over the phase space for different mass ratios. 



2. Examining the set of subchannels in Fig. [I], the QCD corrections can be categorized 
into five classes. Virtual corrections of the internal quark lines, of the 7/Z-quark vertices, 
of the scalar /pseudoscalar-Higgs-quark vertices, and box diagrams interfere with the Born 
amplitude. Gluon radiation off internal and external quark lines adds incoherently to the 
cross sections. In order to disentangle the different subchannels with couplings varying 
over the supersymmetric parameter space, the cross section for scalar Higgs radiation 
must be split into nine parts: Higgs radiation off the quarks (cti,<t 2 ), Higgs-strahlung 
off the Z-boson line ((73,(74), the related interference terms (75, <jq), real/virtual A decay 
(77), and the interference with the first two Higgs-radiation and Higgs-strahlung channels 
(0-8,09): 
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where Hi = h, H. Correspondingly, the cross section for pseudoscalar Higgs radiation 
may be written as 
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split into A radiation off the quarks (o J 1 ,a' 2 ), virtual Higgs decays (a' 3 , a' 4 ) and their 
interference term (cr' 5 ), as well as their interference terms with the radiation diagrams 
(a' 6 ,a' 7 ). The Yukawa and gauge couplings are set off explicitly in (0) and (|3]), with the 
normalized couplings gQQHi etc. given in Table [lfj The electric charges are defined as 
Q e = — 1, Q t = +2/3 and Qb = —1/3, and v ei Q and a e> Q are the vector and axial- 
vector charges of the electron and the quark Q, normalized as v = — 2Q s w)/ S 2W 
and a = Izl/s2Wi with J 3 x = ±1/2 being the weak isospin of the left-handed fermions 
[as usual, sfy = 1 — c 2 ^ = sin 2 6 W = 0.23 and s 2 w — sin20w]; o"o = 47ra 2 /3s denotes 
the standard electromagnetic /i-pair cross section, and N c = 3 is the color factor of the 
quarks. The value of the electromagnetic coupling is taken to be a = 1/128. The mass 
of the Z boson is set to M z = 91.187 GeV, and the pole masses of the top and bottom 
quarks are set to M t = 174 GeV |12] andQ Mb = 4.62 GeV fi~3| , respectively. The masses 
of the MSSM Higgs bosons and their couplings are related to tan /3 and the pseudoscalar 
Higgs boson mass Ma- In the relation we use, higher-order corrections up to two loops 
in the effective-potential approach are included JL4|. The SUSY parameters are chosen as 



pi = A t = Ab = and Mq = 1 TeV; this simple choice is sufficient to illustrate the main 
results. 

The production of bbh/H/A final states in (Q) and (|3p can be mediated by resonance 
channels e + e~ — > Ah/H and Zh/H. We describe the resonance structures as Breit- 
Wigner forms by substituting M 2 — > M 2 — iMT in all boson propagators. The decay 
widths of the Higgs bosons are calculated including higher-order corrections, as described 
in Refs. Hl5jl6"| , while the Z width is set to Tz = 2.49 GeV. For tih/H/A production the 



widths can be neglected, since the Higgs masses are taken below the ti threshold. 

The coefficients cij, in the cross sections (fj) and (0) can be decomposed into Born 
contributions of, erf and QCD corrections Si, 8'f 

= (l + , < = <° (l + ■ (4) 

The renormalization scale of the QCD coupling a s , which is evaluated in NLO with five 
active flavors normalized to a s (M|) = 0.119 |12]1 , is chosen at [i 2 R — s, where s = Eq M 
is the center-of-mass (CM) energy squared. By definition, the Yukawa couplings in the 
Born contributions are evaluated at the scale set by the Higgs momentum flow to leading 
logarithmic order. 

The QCD radiative corrections have been calculated in the standard way. The Feyn- 
man diagrams and the amplitudes for the virtual corrections have been generated with 
FeynArts [|17|]. They have been evaluated by applying the standard techniques for one- 
loop calculations, as described in Refs. |T8| , |T9| . Ultraviolet divergences are consistently 
regularized in D = 4 — 2e dimensions, with 75 treated naively since no anomalies are 
involved. The renormalization of the QQ<p vertices is connected to the renormalization 
of the quark masses, which, in the case of the top quark, is defined on shell (see e.g. 



2 The couplings <7qq</, (gzzc/>, gzA<p) include the normalization coefficients Mq/v {Mz/v), respectively. 
3 This value for the perturbative pole mass of the bottom quark corresponds in NLO to an MS mass 
m b (m b ) = 4.28 GeV. 
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Refs. |19| , |20| for details). Large logarithms are mapped into the running mass ^(Qniggs) 
for the Yukawa couplings of the b quark, with Qm ggs denoting the squared momentum flow 
through the corresponding Higgs-boson line. The algebraic part of the virtual corrections 
has also been checked by using FeynCalc |HJ]. The infrared divergences encountered in 
the virtual corrections and in the cross section for real gluon emission, are treated in two 
different ways. Both calculations follow subtraction procedures, one using dimensional 
regularization and one introducing an infinitesimal gluon mass [22]. The results obtained 
by the two different procedures are in mutual numerical agreement after adding the con- 
tributions from virtual gluon exchange and real gluon emission. A second, completely 
independent calculation of the QCD corrections to the total cross section was based on 
the evaluation of all relevant cut diagrams of the photon and Z-boson self-energies in 
two-loop order, generalizing the method applied to ti(g) intermediate states in Ref. f23fl . 
The results of the two approaches are in numerical agreement. The final expressions for 
the integrated cross sections of, of and for the QCD corrections 6i, 8[ are, however, too 
lengthy to be recorded in this letter^. 



3. The general characteristics of the reactions ([!]) will be discussed for two examples in 
detail, tan/5 = 3 and 30. The numerical analysis has been performed for the case of no 
mixing in the stop and sbottom sectors of the MSSM. Since the top Yukawa couplings are 
suppressed for large values of tan j3, sizeable cross sections for ti(f> production are expected 
only for small and moderate values of tan /3 and in the decoupling regime at large values 
of tan p, where the h couplings approach the SM values. The opposite is realized for bbcf) 
production outside the resonance regions due to the strong enhancement of the b Yukawa 
couplings for large values of tan (3. 

a.) The QCD corrections to the top final states tt(f) can be interpreted easily in two 
kinematical areas. Whenever the invariant mass of the ti pair is close to threshold, the 
gluonic Sommerfeld rescattering-correction is positive and becomes large. In the threshold 
region the K factor approaches the asymptotic form 



with the maximal quark velocity (3 t = y {y/s — M^) 2 — AM 2 /2M t in the (ti) rest frame. 

For high energies, on the other hand, the QCD corrections are of order a s /ir. In the 
energy region s ^> AM 2 ^> Mfj but log s/M 2 ^> 0(1), which is relevant for the present 
analysis, the QCD corrections can qualitatively be traced back to vertex corrections and 
infrared gluon radiation. Since scalar Higgs bosons are radiated off top quarks preferen- 
tially with small energy [x — E$/E t — > 0], as is evident from the leading (universal) part 
of the fragmentation function 



f(t^th/H;x) = ^- 



. 1 — x . s 
4 + x log —2 

x Mf 



(6) 



4 NLO Fortran programs calculating the integrated cross sections and differential distributions for all 
the processes, are available on the web [p4|. 
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the QCD correction of the scalar Yukawa vertex, regularized by soft gluon radiation, 
approaches the value || 



a V+IR _ ^ a s 
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The scalar Yukawa vertex is therefore reduced by four units in a s /7r which are compen- 
sated only partly by one unit due to the increase of the tt production probability, leading 
in total to 

irSt-1-3— + ... f T<j> = h,H. (8) 

7T 

The ellipsis accounts for hard Higgs and gluon radiation (of order +a s /n). Thus, the QCD 
corrections are expected negative for scalar Higgs particles in the high energy continuum. 

By contrast, the corresponding fragmentation function for the pseudoscalar Higgs 
boson 

M~ = P) 

is hard so that the average of the vertex and IR gluon corrections over the Higgs spectrum 
amounts to 

2 — x 
1 H log(l — x) 



A V+IR , 4 "s 



(10) 

X ' J / 2 7T V ^ 

Adding to this correction the increase of the ti production probability of one unit, the K 
factor is very close to unity 

*£i-i-^ + ... (ii) 

After hard gluon bremsstrahlung is taken into account (symbolized by the ellipsis), the 
overall QCD corrections for the pseudoscalar Higgs boson are therefore expected slightly 
positive. [For ultra-high energies, i.e. logs/M 4 2 3> 1, hard gluon bremsstrahlung becomes 
important. Similarly to the leading terms in the fragmentation functions eqs. (|6|) and (|9]), 
the QCD corrections for scalar and pseudoscalar Higgs bosons approach each other as a 
result of chiral symmetry restoration in asymptotia; this has been verified in a numerical 
calculation.] 

Similar estimates can be applied to bottom final states which in general are dominated 
by resonance decays. After absorbing the large logarithms log^^/M 2 ) into the Yukawa 
couplings, the non-leading effects are positive: 



+ for {Higgs, Z}^ bb . (12) 

7T 



Also close to the thresholds and in the high-energy limit the QCD corrections remain 
positive after mapping the large (negative) corrections into the running Yukawa couplings. 
Since different channels are activated at the same time, only a qualitative estimate can 
be given in the continuum regime, 

<! = l + c- withe = 0(1), (13) 

7T 
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while details must be left to the numerical analysis. 

b.) The results are exemplified in Fig. |2| for tan/3 = 3 and 30 and for the collider 
energy E C m = 500 GeV. If required by the size of the cross section, which should not fall 
below ~ 10~ 2 fb in order to be accessible experimentally, we switched to Eqm = 1 TeV. 
The Born terms are shown by the dotted curves, while the final results for the cross 
sections, including QCD corrections, are given by the full curves. 

For E C m = 500 GeV the QCD corrections to Higgs-boson production in association 
with tt pairs increase the scalar Higgs-pro duct ion cross sections significantly, as can be 
inferred from Fig. Close to threshold the numerical results clearly exhibit the strong 
increase of the cross sections due to the Coulomb singularity @. Moreover, for tan/3 = 30 
the cross sections are strongly suppressed except for the regions where the light (heavy) 
scalar Higgs mass is close to its upper (lower) bound. For tan/3 = 3 the cross section 
amounts to about 1 fb, which leads to a significant number of events at the TESLA 
collider, being designed to reach integrated luminosities of about J C ~ 1 ab~ x in three 
years of operation. 

For CM energies of 1 TeV the QCD corrections to scalar Higgs production in associa- 
tion with tt pairs are of moderate size. They decrease the cross sections by about 3-5%. 
This is in accordance with the asymptotic form of the K factor |7J . For pseudoscalar Higgs 
production, the size of the QCD corrections is slightly positive at 1 TeV, in agreement 
with the qualitative discussion above. 

The cross sections for Higgs-boson production associated with bb pairs are significantly 
larger due to the resonance contributions from on-shell Z and Higgs-boson decays into 
bb pairs. The QCD corrections increase these cross sections by about 5-25%. The drop 
in the bbA production cross section for tan (3 = 3 at Ma ~ 140 GeV and 175 GeV can be 
attributed to the crossing of the thresholds for resonant H — > WW and H — > hh decays, 
respectively, in HA final states. 

Without cuts, the intermediate resonance decays Z, h,H,A —>■ bb dominate all bbcf) 
production processes, whenever they are kinematically allowed, and it will be difficult 
to extract the bottom Yukawa couplings experimentally in regions where resonant Higgs 
decays to bb pairs are dominant. This is the case at large values of tan/3 for all neutral 
Higgs particles and at small values of tan/3 for Higgs masses below the WW threshold. 
In these cases the branching ratios, which determine the size of the bbcf) cross sections, 
will be nearly independent of the bottom Yukawa couplings. Resonance decays R — > bb 
pairs in the bbh/H/A final states can however be eliminated by cutting out the resonance 
energy of the final-state Higgs boson, — (s + M| — Mj { )/2y/s. After subtracting 

these resonance parts, the non-resonant contributions are suppressed by about one to 
three orders of magnitude. The resonances have been removed from the cross sections in 
the examples of Fig. ||] by subtracting the two-particle cross sections in the Breit-Wigner 
bands Mr ± A of the energy E^^ with the resolution A = 5 GeV. This theoretical 
definition is used for the sake of simplicity; wider cuts may be required in experimental 
analyses. Peaks and dips in the cross sections are the result of overlapping Breit-Wigner 
bands. For the scalar Higgs particles they arise from overlapping Z and A boson bands; 
in pseudoscalar Higgs production the light and heavy scalar resonance bands overlap for 
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100 GeV < Ma ^ 120 GeV for tan/5 = 30. The dips occur whenever the two resonance 
bands touch each other, while the peaks between the dips occur when the resonance 
masses coincide exactly. As shown in Fig. |3], the QCD-corrected cross sections are still 
close to lib or slightly below, except for heavy masses at small values of tan/5. Thus, 
ensembles of order 10 3 events can be collected at a high-luminosity collider. 

4. Measurements of Yukawa couplings in supersymmetric Higgs radiation off heavy quarks 
at e + e~ linear colliders are difficult. This is a result of the large number of subchannels 
contributing to the QQh/H/A final states in supersymmetric theories in general, and 
the contamination by two-Higgs final states in particular. Nevertheless, the continuum 
cross sections appear large enough to approach this experimental problem as shown in 
the present analysis. Even though experimental simulations are beyond the scope of this 
note, it may be concluded from earlier SM analyses that the method will work at least 
in part of the MSSM parameter space, thus providing us with the absolute size of the 
quark-Higgs Yukawa couplings in the minimal supersymmetric theory. 
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Figure 2: Production cross sections for the MSSM Higgs bosons h, H and A in association 
with heavy t, b quark pairs; Born approximation: dashed, QCD corrections included: full 
curves. The rapid drops in the bbA cross section at tan (3 = 3 are due to the kinematical 
opening of the resonant H — > WW, hh decays, which reduce the branching ratio of the 
resonant H — > bb decay. 
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M h/H [GeV] M A [GeV] 

Figure 3: Continuum production of MSSM Higgs bosons h, H, A in association with a bb 
pair after subtraction of resonance decays to bb pairs in the Breit- Wigner bands Mr ± 5 
GeV. 
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